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Abstract. The uncertainties in Parton Distribution Functions (PDFs) are the dominant source of the 
systematic uncertainty in precision measurements of electroweak parameters at hadron colliders (e.g. 

0eff{Mz), Ow = 1 — M^/M% and the mass of the W boson). We show that measurements 
of the forward-backward charge asymmetry [ApB^M^y)) of Drell-Yan dilepton events produced at hadron 
colliders provide a new powerful tool to reduce the PDF uncertainties in these measurements. 

PACS. 12.5.-y electroweak 12.38.-t Quantum chromodynamics 


1 Introduction 

Precision measurements in hadron colliders are limited by 
our knowledge of Parton Distribution Functions (PDFs). 
In general, PDF fits by various groups including cteq [T], 
mmht[2], nnpdf [3l|4], hera[5], and abm[6] are extracted 
from fixed target experiments and various cross sections 
measurements at colliders. The fixed target experiments 
include electron, muon, neutrino, and Drell-Yan experi¬ 
ments. The collider experiments include ep(HERA), pp 
(Tevatron) andpp(LHC). 

Some of the fixed target measurements are on nuclear 
targets resulting in additional uncertainties from modeling 
of nuclear effects. Some of the fixed target measurements 
are also at low momentum transfers where the contribu¬ 
tions of non-perturbative and higher twist effects may be 
significant. These issues are absent in collider cross section 
data. Therefore, recent PDF fits have placed a greater em¬ 
phasis on collider cross section data. 


1.1 Measurements of electroweak parameters at 
hadron colliders 

Within the standard model, measurements of the mass 
of the Z boson and top quark, in combination with the 
mass of the Higgs boson, can be used to predict the mass 
of the W boson (Mw). At present, the average of the 
all direct measurements of Mw (80385T15 MeV) is about 
1.5 standard deviation higher [7] than the prediction of the 
standard model. Predictions of supersymmetric models for 
Mw are also higher than the predictions of the standard 
model [13]. Therefore, more precise measurements of the 
mass of Mw are of great interest. 


Alternatively, Mw can also be extracted indirectly from 
measurements of the on-shell electroweak mixing angle 
sin^ Ow by the relation sin^ Ow = 1 — ^wl^z- 

Measurements of the forward-backward charge asym¬ 
metry in Drell-Yan dilepton events produced at hadron 
colliders (in the region of the Z pole) have been used to 
measure the value of the effective electroweak (EW) mix¬ 
ing angle sin^ jOlllOillfHT^ . In addition, by incor¬ 

porating electroweak radiative corrections in the analysis 
the CDE collaboration has also measured the on-shell EW 
mixing angle sin^ Ow |9illQj • 

An uncertainty of ±0.00030 in the measurement of 
sin^ Ow is equivalent to an indirect measurement of Mw 
to a precision of ±15 MeV. However, the PDE uncer¬ 
tainty quoted in the most recent measurement of sin^ O^f f 
by the ATLAS collaboration [T2] at the LHC is ±0.00090. 
Therefore, a significant reduction in the PDE uncertainty 
is needed. In this communication, we show how Afb data 
also provide a new powerful tool to reduce PDE uncer¬ 
tainties in the measurements of electroweak parameters 
in hadron colliders 

The constraints provided by Afb measurements in 
combination with constraints from the W charge asym¬ 
metry (Aw) can be used to reduce the PDE uncertainty 
in the extracted value of sin^ Ow and sin^ 0^fJ^^{Mz) from 
Afb data. The Afb constraints on PDEs can also be used 
to reduce the PDE uncertainty in other precision measure¬ 
ments with Z and W bosons such as the measurement of 

Asymmetries such diS Afb and Aw are ideal in pro¬ 
viding additional constraints because asymmetries are less 
sensitive to the choice of QCD scale and QCD higher order 
terms. In addition, there are new techniques that can be 
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used p^fTS] to greatly reduce the experimental systematic 
uncertainty in asymmetry measurements. 


2 qq annihilations to dileptons 


In leading order (LO) dileptons are primarily produced in 
quark-antiquark annihilation. Here, one parton (quark or 
antiquark) carries momentum xi and another parton car¬ 
ries momentum X 2 . The momentum fractions Xi ^2 carried 
by the partons are related to the mass (M) and rapidity 
(y) of the two leptons as follows: 


M 

xi^2 = 

VS 




(1) 


The angular dependence of the differential cross section 
for qq annihilation to a dilepton pair can be written as 



M (GeV/c2) 

Fig. 1. The contributions of u-type quarks (blue) and d-type 
quarks (red) to Afb{M) at the Tevatron. 


oc (1 + cos^ 0) + AAM) cos 0 (2) 

a(cos 0) 

where 0 is the emission angle of the negatively charged lep¬ 
ton relative to the quark momentum in the dilepton center 
of mass frame, and A 4 (M) is parameter that depend on 
the weak isospin and charge of the incoming quarks. 

The cross sections for forward {(Jf) and backward {(Jb) 
events are given by 

«+ 0(0 

The electroweak interaction introduces an asymmetry 
(a linear dependence on cos6>), which can be expressed as 

Afb{M) = = \Ai{M) (5) 

(Tp ^ CTb O 

The dependence of ApsiM, y) on sin^ has been used 
to measure sin^ at the Tevatron and LHC. 

The systematic uncertainties in the measurement of 
ApBiM^y) can be greatly reduced if ApBiM^y) is ex¬ 
tracted from a measurement of A 4 ^{M^y). This can done 
by an event weighting technique [15] for which there is a 
cancelation of systematic errors that originate from un¬ 
certainties in acceptance and efficiencies. With this tech¬ 
nique, no acceptance or efficiency corrections are needed. 
The extracted values of A/^i^M^y) using the event weight¬ 
ing technique are equal to the Born level A/^i^M^y). This 
technique has been using in the most recent measurements 
at CDF [10] . 


3 Afb at the Tevatron 


For pp collisions, the direction of the quark is predomi¬ 
nately in the proton direction, and the direction of the an¬ 
tiquark is predominately in the antiproton direction. Here, 
most of the cross section originates from the annihilation 
of quarks in the proton with antiquarks in the antiproton. 
Therefore, Afb is measured under the assumption that 
the quarks originate form the proton, and the antiquarks 
originate from the antiproton (first term in eqj^. 

Since q{x) in the proton is equal to q{x) in the antipro¬ 
ton, the dilepton production cross section can be expressed 
as follows: 

^(pp)oc Vi{qi{xi) ■ qi{x2) + Qiixi) ■ qi{x2)} (Q) 

flavor 


Here qi{x) denote the quark distributions {u{x)^ d{x)^ 
s{x) , c(x), b{x)) and qi{x) denotes the antiquark distribu¬ 
tions ('u(x), d(x), s{x), c(x), b{x)) in the nucleon. The pa¬ 
rameters Vi denote the Z/y couplings for each flavor. Here, 
Vi are functions of both the dilepton mass and sin^ . 

The extraction of sin^ from Afb{M) (or A^M)) 
is sensitive to PDFs for two reasons. First, Afb{M) for 
charge 2/3 (u-type) quarks and charge 1/3 (d-type) quarks 
is different. Fig.[^shows the contributions of u-type quarks 
(blue), d-type quarks (red) and the sum of the two con¬ 
tributions (black) to Afb{M) at the Tevatron as given 
by 

j^d—type ^ {dcTjF {dd) b 

(dd) F + {dd) B + {uu) F + {uu) b 
j^u—type ^ {uV^F {uV) B 

{dd)F + {dd)B + {uu)f + {uu)b 

The measured asymmetry is sensitive to the fraction of 
down quarks in the proton because the asymmetries for 
up and down quarks are different. The sensitivity is pro¬ 
portional to 


D 


T ev 
AFB 


{d)9S 


d{xi) d{x 2 ) 
u{xi) u{x 2 ) 




(7) 
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In addition, there is a small fraction of events for which 
the annihilation is between sea antiquarks in the proton 
with a sea quarks in the antiproton (second term in eq.[^. 
The forward-backward asymmetry Afb{M) of the second 
term in equationl^is opposite to the Afb{M) of the larger 
first term. This also results in a dilution of the 

measured asymmetry. 


D 


Tev ( 
AFB 


(q) 


T.flavor Vqq{Xl) ■ q{X 2 ) 


U{xi)u{x2) 


( 8 ) 


The antiquark dilution is primarily from u type anti¬ 
quarks. For proton-antiproton collisions, most of the cross 
section is near y=0 {xi ~ X 2 ). Therefore, the PDF un¬ 
certainty in the extraction of sin^ from Afb{M) (or 
A 4 {M)) at the Tevatron depends primarily on how well we 
can constrain the following contributions to the dilution 
at xi = 


PDFs. Each of the PDF replicas has equal probability of 
being correct. The central value of any observable is the 
average of the values Si= {sin^ Ow)i extracted with each 
one of the N PDF replicas. The PDF uncertainty {=(Tpdf) 
is the RMS of the values extracted using all N replicas. 


1 ^ 

i=l 

(12) 

^ _jFlAsi-{sW 

pdf Y _i 

(13) 


and the uncertainty in the estimate of the PDF uncer- 
tainty is Aapdf = 

The two methods provide equivalent information. For 
any given a set of Hessian eigenvector PDFs there is a 
prescription to generate j^llSillQj an arbitrary number of 
PDF replicas. 


Dl^Md) « [^(^ 1 )]" (9) 

Dl7B{q) « [^(^ 1 )]" ( 10 ) 


3.1 W charge asymmetry at the Tevatron 

The W~ /W^ ratio at the Tevatron can be written as 

/ VTev ^ ^(^1) ^^(^2) + <g((^l) C{X2) _ d d 

IU+ u{xi) d{x 2 ) c{xi) s{x 2 ) u u 

( 11 ) 

Precise measurements of the W asymmetry provide 
information on the d/u ratio at the Tevatron. These mea¬ 
surements are important to constrain the PDF uncertain¬ 
ties for the direct measurement of the W mass. However, 
at the Tevatron these measurements do not provide infor¬ 
mation relevant to the measurement of sin^ Oef / for two 
reasons. First, there is no information at y=0 {xi ~ X 2 ) 
since here the W charge asymmetry at the Tevatron is 
zero. Secondly, at the Tevatron, the W charge asymme¬ 
try does not provide information on the absolute level of 
^ {x ). The W charge asymmetry at the Tevatron provides 
information only on the slope of ^{x) as a function of x. 


3.2 PDF uncertainties: Hessian and Replica PDFs 

All PDF groups provide a default (central) PDF set. There 
are two methods that are used for the determination of 
PDF uncertainties. The first method is to provide a set 
of eigenvector error PDFs (Hessian method). The PDF 
uncertainties in a measurement are determined by repeat¬ 
ing the analysis for all of the error PDF sets, and adding 
in quadrature the difference in the results obtained with 
the error PDFs and the results obtained with the default 
PDF. 

The second method (which is referred to as replica 
PDFs) is to provide a set of N (e.g. 100 or 1000) replica 


3.3 Reducing PDF uncertainties with new data 

The advantage of the PDF replica method is that con¬ 
straints from new data can easily be incorporated in any 
analysis by applying different weights for each replica. 

Replicas for which the theory predictions are in agree¬ 
ment with the new data are given higher weights, and 
replicas for which the predictions are in poor agreement 
are given lower weights. The weights are derived from the 
values of the comparison between the new data and 
theory prediction each of the PDF replicas. 

The central value of any observable is the weighted 
average of the values extracted using each one of the N 
PDF replicas. The PDF uncertainty is the weighted RMS 
(root mean square) of the values extracted each of the N 
replicas. 

The procedure of including constraints from new data 
was initially proposed by Giele and Keller [20]. They pro¬ 
posed that each of the N PDF replicas be weighted as 
follows: 


Wi 

e 2 Vi 

Ei=i e 

(14) 


N 


{s) 

= y^^WjSj 

(15) 


i=l 


^pdf 

A:l,w4si-{s)r 

V 1 - We// 

(16) 


The weights reduce the effective number of replicas 
from N to A^e// where 




1 


eff 


E N 2 

i=i < 


(17) 


and the uncertainty in the estimate of the PDF uncer¬ 
tainty is Aapdf — - 
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Fig. 2. Afb versus dilepton mass at the Tevatron for 
sin^ ^h^= 0.2244 and the default NNPDF 3.0 (nnlo) PDF 
(261000). The band corresponds to ten NNPDF replicas. 

More recent discussions of the method can be found in 
references p!8l[T^[^[l22[|23j . In the sections that follow we 
show how the mass and rapidity dependence oi Afb can 
be used to both provide additional constraints and reduce 
the PDF uncertainty in measurements of sin^ 


3.4 Number of replicas needed 

Typically between 100 and 1000 PDF replicas are used. 
A large number of replicas is only needed if the new data 
that is being incorporated is so precise that the number of 
effective replicas drops below 10. This only happens if the 
statistical errors of the new data are much smaller than 
the PDF uncertainties. 

For the electroweak measurements that are discussed 
in this paper the statistical errors which are achievable 
in the next few years are typically within a factor of 2- 
3 of the PDF uncertainties. Therefore, 100 replicas are 
typically sufficient. 

3.5 Mass dependance of Afb{M) as a function of 
sin^ Ow and PDFs at the Tevatron 

The sensitivity of the mass dependence of Afb{M) on 
sin^ Ow and PDFs is different. In the region of the Z pole, 
Afb{M) is sensitive to the vector couplings, which depend 
on sin^ Ow- At higher and lower mass Afb{M) is sensitive 
to the axial coupling and therefore insensitive to value of 
sin^ 0w‘ 

In contrast, the magnitude of the dilution of Afb{M) 
depends on the PDFs. The sensitivity to PDFs is largest 
in regions where Afb{M) is large (i.e. away from the Z 
pole). 

Fig. shows Afb{M) as a function dilepton mass at 
the Tevatron for sin^ 6 >vf= 0.2244. The band corresponds to 
the predicted values of Afb{M) for the default NNPDF 3.0 
(nnlo) pdf (261000), and ten NNPDF 3.0 (nnlo) replicas. 
Afb{M) is shown for y^=1.96 TeV and dilepton rapid¬ 
ity less 1.7, which corresponds to a typical acceptance for 
Tevatron experiments (CDF or DO). 


Fig.j^a) shows the sensitivity of Afb{M) at the Teva¬ 
tron to PDFs. The lines are the difference between Afb{M) 
for 10 NNPDF 3.0 (nnlo) replicas and Afb{M) calculated 
for the central dekult NNPDF 3.0 (nnlo) (261000). Here 
sin^ Ow is fixed at a value of 0.2244. The difference origi¬ 
nates from the differences in ^{x) and the antiquark frac¬ 
tions for the different PDF replicas. 

Fig.j^b) shows the sensitivity of Afb{M) at the Teva¬ 
tron to sin^ Ow • The lines are the difference between the 
calculated Afb{M) for sin^ Ow values ranging from 0.2220 
(show at the top in red) to 0.2265 (shown in the bottom 
in blue) and Afb{M) for sin^ ^vf= 0.2244. Here Afb{M) 
is calculated with the default NNPDF 3.0 (nnlo) (261000). 




Fig. 3. Tevatron: (a) The difference between Afb{M) for 
10 NNPDF 3.0 (nnlo) replicas and Afb{M) calculated for the 
default NNPDF 3.0 (nnlo) (261000). Much of the difference 
originates form the different dilution factors for each of the 
NNPDF replicas. Here sin^ Ow is hxed at a value of 0.2244. (b) 
The difference between Afb{M) for different values of sin^ Ow 
ranging from 0.2220 (shown at the top in red) to 0.2265 (shown 
on the bottom in blue), and Afb{M) for sin^ ^vf=0.2244. Here 
Afb{M) is calculated with the default nnpdf 3.0 (nnlo). 

As shown in Fig.|^a) there is a large difference in the 
Afb{M) predictions for PDF sets with different ^{x) and 
antiquark fractions |(x) in regions where Afb{M) is large 
and positive (M>100 GeV). The changes in Afb{M) in 
regions where Afb{M) is large and negative (M<80 GeV) 
are in the opposite direction. 

In contrast, as shown in Fig. [^b), different values of 
sin^ Ow change Afb{M) primarily in the region near the 
Z pole. However, here the change is in the same direc¬ 
tion above and below the Z pole. Therefore, if we extract 
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sin^ Ow from Afb{M) data with different PDFs, PDFs 
with poor values of are less likely to be correct. 


3.6 MC studies of dilepton production at Tevatron 

The 10 fb“^ Run II e+e“ data sample at CDF corresponds 
to about 500K events. A similar sample was collected by 
the DO experiment [TT]. The acceptance of the Tevatron 
experiments limits the sample to events with dilepton ra- 
pidity I 2 /I <1.7. 

We simulate Afb{M) measurements corresponding a 
10 fb“^ statistical sample at the Tevatron with three dif¬ 
ferent input assumptions for Afb^ Iri all cases we use 
sin^ =0-2244 and calculate Afb in 15 bins for dilep¬ 
ton mass spanning the range from M=50 GeV to M=150 
GeV. We generate pseudo data for three input assump¬ 
tions. For each input assumption we generate a set of 1600 
pseudo-experiments. 

— The input assumption for the first set of 1600 pseudo 
experiments is that Afb{M) is equal to the predic¬ 
tions of a Tree-level calculation (including EBA EW 
radiative corrections [OllTO] i calculated with the default 
NNPDF 3.0 (nnlo) PDE set. 

— The input assumption for the second set of 1600 pseudo 
experiments is that Afb{M) is equal to the predic¬ 
tions of a Tree-level calculation (including EBA EW 
radiative corrections [QllTO ] ) calculated with the default 
NNPDF 2.3 (nnlo) PDE set. 

— The input assumption for the third set of 1600 pseudo 
experiments is that Afb{M) is equal to the predictions 
of RESBOS [16] (modified to include EBA EW radiative 
corrections [OlfTO] i calculated with the cteq 6.6 PDE 
set. 



Fig. 4. An example of the extraction of sin^ Ow from Afb{M) 
data at the Tevatron. Here, XAfb is plotted for different values 
of sin^ Ow The extracted value of sin^ Ow is the value with 
the minimum XAfb the statistical error corresponds to a 
change of by ±1. 


Table 1. Values of sin^ Ow with statistical errors and PDF un¬ 
certainties expected at the Tevatron for a 10 fb“^ sample for a 
CDF like detector. The PDF uncertainty for a standard anal¬ 
ysis is compared to the PDF uncertainty for an analysis with 
XAfb weighting. The default NNPDF 3.0 (nnlo) is used to gener¬ 
ate the pseudo data in the first column and the default NNPDF 
2.3 (nnlo) is used to generate the pseudo data in the second 
column. All pseudo data are generated with sin^ =0.22420. 


CDF-like detector 
Pseudo-Experiment 
Tevatron 10 fb~^ 

500A reconstructed 
e^e“ events 

Input 
Tree-level 
Default 
nnpdf 3.0 
(261000) 
(nnlo) 

Input 
Tree-level 
Default 
nnpdf 2.3 
(261000) 
(nnlo) 

sin^ Ow input 

0.22420 

0.22420 

statistical error 

A sin^ Ow 

CTIO PDF error 

±0.00042 

±0.00026 

±0.00042 

±0.00026 

Number of analysis replicas 
nnpdf replica set 
Templates 

100 

nnpdf 3.0 
Tree-level 

100 

nnpdf 2.3 
Tree-level 

Average method 
extracted sin^ Ow 

PDF error RMS 
(uncertainty in PDF error) 

Neff = 100 

0.22420 
±0.00027 
( 0.00002) 

Neff = 100 

0.22420 
±0.00028 
( 0.00002) 

XAfb weighting 
extracted sin^ Ow 

PDF error weighted 
(uncertainty in PDF error) 

Neff — 88 

0.22420 
±0.00020 
( 0.00002) 

Neff =85 

0.22420 
±0.00022 
( 0.00002) 


3.6.1 Tevatron pseudo data: default nnpdf 3.0 (nnlo) 
and default nnpdf 2.3 (nnlo) 

For the first set of 1600 pseudo experiments the default 
NNPDF 3.0 (nnlo) is used to generate pseudo data. The 
simulated values of Afb{^) for each experiment are com¬ 
pared to Afb{M) templates generated at Tree-level for a 
range of values of sin^ Ow for each of the 100 nnpdf 3.0 
(nnlo) pdf replicas. For each replica we extract the best 
fit value of sin^ Ow^ the corresponding statistical error and 
the fit XAfb’ There are about 500K dimuon events in each 
Tevatron pseudo-experiment, which results in a statistical 
error in sin^ Ow of ±0.00042, 

An example of the extraction of sin^ Ow from Afb{M) 
data at the Tevatron is shown in Fig0 Here, Xa/i, is plot- 
ted for different values of sin^ Ow The extracted value of 
sin^ Ow is the value with the minimum x^fb the sta¬ 
tistical error corresponds to a change of XAfb hy ±1- 

For the second set of 1600 pseudo experiments the de¬ 
fault NNPDF 2.3 (nnlo) is used to generate pseudo data 
and the extraction of sin^ Ow is done using 100 nnpdf 2.3 
(nnlo) pdf replicas.. 

For each set, the extracted value of sin^ Ow and the 
PDF uncertainty are done in two ways. 
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1. the standard average and RMS of the sin^ Ow values for 
the 100 PDF replicas. 

2. the weighted average and weighted RMS of the 

sin^ Ow values for the 100 PDF replicas. 


^NNPDF3.0 Default pseudo data (stat error) 

24 

22 

c 

CM E 20 

18 

16 

14 

0.2: 

;100 NNPDF3.0 Replicas 

® ■ NNPDFaveraaefPDFerror^ 


• 

• 

• 

• 

• • 

*• J • • n-‘* 

• ••• 

232 0.2237 0.2242 0.2247 0.2252 

sin^e^ 


^NNPDF2.3 Default pseudo data (stat error) 


100 NNPDF2.3 Replicas 

24 



22 

^ weiehtine (PDF error! 


c 

• 

1 X 

• 


• 

18 

• - . • • • 

*.•••• , • 



16 


. 

14 

1 1 1 1 1 1 1 1 1 1 

0.2232 0.2234 0.2236 0.2238 0.224 0.2242 0.2244 0.2246 0.2248 0.225 0.2252 


sin^Ow , 


Fig. 5. Tevatron: A graphical illustration of the analysis 
of one typical pseudo experiment. Shown is a scatter plot of 
Ow and Xa/ 6 values for 100 PDF replicas, (a) For pseudo 
experiment generated with the default NNPDF 3.0 (nnlo) and 
sin^ =0.22420 at Tree-level, (b) or pseudo experiment gener¬ 
ated with the default NNPDF 2.3 (nnlo) and sin^ ^vf=0. 22420 
at Tree-level. Also shown on the plot is the input value of 
sin^ Ow with the average statistical error of one pseudo experi¬ 
ment. In addition, we show the average of the extracted values 
sin^ Ow and average PDF uncertainty for both the standard 
analysis, and the Xa/ 6 weighted analysis. 


For each of the 100 NNPDF 3.0 (nnlo) (or NNPDF 2.3 
(nnlo)) replicas we calculate the average of the 1600 ex¬ 
tracted values of siu^ Ow^ the average of the 1600 PDF uu- 
certaiuties, aud the average 1600 statistical errors. These 
average quautities have suiall fiuctuatiou aud represeut 
the result of oue pseudo experiuieut ou average. The av¬ 
erage of the 1600 PDF uucertaiuties is au estiuiate of the 
typical uucertaiuty for oue iudividual pseudo experiuieut. 
lu order to test for possible bias iu the uiethod, the av¬ 


erage of the 1600 extracted values of siu^ Ow is couipared 
the 0.22420, which is the value used iu the geueratiou. 

As expected iu both aualyses the average extracted 
value of siu^ Ow is the sauie as the value with which the 
pseudo data has beeu geuerated (0.2242), as showu iu Ta¬ 
ble With the weightiug uiethod the PDF uucer¬ 
taiuty iu the extracted value of siu^ Ow is reduced froiu 
±0.00027 to ±0.00020. This illustrates that although the 
statistical error iu siu^ Ow of ±0.00042 is souiewhat larger 
thau the PDF uucertaiuty of ±0.00027, the Afb data at 
higher aud lower luass has sufficieut precisiou to coustraiu 
the PDFs which yields a 25% reductiou iu the PDF uu¬ 
certaiuty. 

A graphical illustratiou of the uiethod is showu iu iu 
Fig.[^(a) aud Fig. [^(b). For each PDF replica, we calcu¬ 
late the average of the extracted values of siu^ Ow and the 
average XAfb i^de fits for the 1600 pseudo experiiueuts. 
Fig. (a) aud Fig. (b) show the scatter plot of the av¬ 
erage of the extracted values of siu^ Ow and the average 
XAfb i^de 100 PDF replicas. 

Also showu ou the plot is the iuput value of siu^ Ow 
with the average statistical error of oue pseudo experi¬ 
uieut. Iu additiou, we show the average of the extracted 
values siu^ Ow aud average PDF uucertaiuty for both the 
staudard aualysis, aud the XAfb weighted aualysis. 


Table 2. Tevatron Pseudo data generated with RESBOS and 
CTEQ 6.6 PDFs, for a CDF like detector. Here, we compare 
sin^ Ow values with statistical errors and PDF uncertainties ex¬ 
tracted with NNPDF 3.0 (nnlo) pdf replicas and with NNPDF 
2.3 (nnlo) pdf replicas. The values of sin^ Ow extracted with 
NNPDF 3.0 (nnlo) pdf and nnpdf 2.3 (nnlo) are different 
(for details see text). 


Pseudo-Experiment 
Tevatron 10 fb“^ 

500K reconstructed 
e^e“ events 

RESBOS 

CTEQ 6.6 

RESBOS 

CTEQ 6.6 

sin^ Ow input 

0.22420 

0.22420 

statistical error 

±0.00042 

±0.00042 

A sin^ Ow 

CTIO PDF error 

±0.00026 

± 0.00026 

Number of analysis replicas 

100 

100 

NNPDF replica set 

NNPDE 3.0 

NNPDF 2.3 

Templates 

Tree-level 

Tree-level 

Average method 

Neff = 100 

Neff = 100 

extracted sin^ Ow 

0.22425 

0.22469 

bias 

±0.00005 

±0.00049 

PDF error RMS 

±0.00027 

±0.00027 

(uncertainty in PDF error) 

( 0.00002) 

( 0.00002) 

Afb XAfb weighting 

Neff = 88 

Neff = 63 

extracted sin^ Ow 

0.22425 

0.22452 

bias 

±0.00005 

±0.00032 

PDF error weighted 

±0.00020 

±0.00021 

(uncertainty in PDF error) 

( 0.00002) 

( 0.00002) 
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Fig. 6. Analysis of a Tevatron pseudo-experiment. The 
pseudo data are generated by RESBOS with CTEQ 6.6 PDF 
and sin^ =0.22420. This figure illustrates that with the x\fb 
weighting method we can determine that pseudo data gener¬ 
ated with CTEQ 6.6 PDFs are not consistent with the NNPDF 2.3 
(nnlo) set. (a) Analysis with 100 nnpdf 3.0 (nnlo) replicas, 
(b) Analysis with 100 NNPDF 2.3 (nnlo) replicas. The distribu¬ 
tion of x\fb values versus sin^ Ow provides a powerful tool to 
discriminate against PDF sets which are incompatible with the 
data. The PDF sets which are compatible with the data should 
have a symmetric distribution of Xa fb values versus sin^ Ow • 


3.6.2 Pseudo data: resbos with cteq 6.6 PDF set 


In the analysis of the RESBOS /CTEQ 6.6 pseudo data 
with NNPDE 3.0 (nnlo) replica templates we find that the 
PDF uncertainty in the extracted value of sin^ Ow when 
we use the standard average is ±0.00027. The PDF un¬ 
certainty is reduced to ±0.00020 when the XAfb weighting 
method is used, as shown in Table and Fig. The ef¬ 
fective number of replicas is reduced from 100 to 88. The 
average value is sin^ Ow = 0.22425 for both the standard 
analysis and the XAfb weighting analysis. The very small 
difference (±0.00005) from the input value of sin^ Ow = 
0.22420 is attributed to the difference between the RESBOS 
pseudo data which is generated at NLO and the templates 
which were done at LO Tree-level. 

In contrast, the standard analysis with the nnpdf 2.3 
(nnlo) replica templates yields a value which is biased 
by ±0.00049±0.00001. This is larger than the PDF un¬ 
certainty of ±0.00027. This bias indicates that the nnpdf 
2.3 (nnlo) set is not fully consistent with the cteq 6.6 
PDF for the Bjorken x region for the production of Z 
bosons at the Tevatron. When the XAfb weighting tech¬ 
nique is used instead, the bias is partially reduced from 
±0.00049±0.00001 to ±0.00032±0.00001, and the effec¬ 
tive number of PDFs is reduced from 100 to 63. The re¬ 
duced bias is expected because XAfb weighting assigns 
small weights to a fraction of nnpdf 2.3 (nnlo) PDF 
replicas which are incompatible with the cteq 6.6. pseudo 
data. 

As shown in Fig. the distribution of XAfb values 
versus sin^ Ow provides a powerful tool to discriminate 
against PDF sets which are incompatible with each other 
or with the data. Our study indicates that cteq 6.6 PDFs 
are inconsistent with the nnpdf 2.3 (nnlo) set, but are 
consistent with the nnpdf 3.0 (nnlo) set. One of the dif¬ 
ference between nnpdf 3.0 and nnpdf 2.3 is that nnpdf 
2.3 used W asymmetry data which is now known to be 
incorrect. 


4 Production of dilepton events at the LHC 


At the LHC, dileptons are produced by annihilation of 
quarks in one proton with antiquarks in the other proton. 


da 

dMdy 


(pp) oc 


T] + qi{xi)qi{x2)} (18) 

flavor 


We perform two analyses of the third set of 1600 pseudo 
experiments (cteq 6.6 pseudo data). In one analysis the 
simulated values of Afb{M) for each experiment are com¬ 
pared to templates calculated at Tree-level for each of the 
100 NNPDF 3.0 (nnlo) PDF replicas. In the other analysis 
the simulated values of Afb{M) for each experiment are 
compared to templates calculated at Tree-level for each of 
the 100 NNPDF 2.3 (nnlo) PDF replicas. In each of the 
two analyses, sin^ Ow is extracted using both the standard 
average and RMS , and also the x\fb weighted average and 
RMS of the 100 PDF replicas. The results are summarized 
in Table [21 


Because on average, quarks carry more momentum than 
antiquarks, the quark direction is assumed to be the di¬ 
rection of motion of the dilepton pair. This is more likely 
to be true for dileptons produced at high rapidity. At the 
LHC the asymmetry from the first term of equation [T8| is 
diluted by the asymmetry of the second term (which is in 
the opposite direction). Equation 18 shows that for y=0 
{ xi = X 2 ) the asymmetries for the two terms cancel each 
other. 

An estimate of the dilution of Afb{M) can be ob¬ 
tained from the probability to misidentify the direction of 
the quark f{M, y). For pp collisions f{M, y) is the fraction 
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Fig. 7. LHC: Top panel: Afb at the LHC at v^=8 TeV 
for six rapidity bins (iY=0 to 5) with average \y\ values of 
0.2, 0.6, 1.0, 1.4, 1.8 and 2.2. For each rapidity bin there are 
twelve mass bins discussed in the text (iMass=0 to 11). The 
horizontal scale for each of the six plots is the dimuon invari¬ 
ant mass for each rapidity bin expressed as 12 x zY + iMass. 
Bottom panel: The green bands span the difference between 
Afb{M) calculated for the 100 nnpdf 3.0 (nlo) replicas and 
Afb{M) calculated for the central default NNPDF 3.0 (nlo) for 
the six dimuon rapidity bins. The blue lines are the differences 
between Afb{M) calculated with different values of sin^ Oeff 
(0.23120 ±0.00040, ±0.00080 and ±0.00120). and the values 
calculated with nominal sin^ ^e//=0.23120. For all of the blue 
lines, Afb{M) is calculated with the central default NNPDF 
3.0 (nlo). The calculations are done with the POWHEG MC 
generator. 



of events for which the antiquark carries more momentum 
than the quark. 




Y.flavor MQi(xl) qi(X2)} 

YflavorM<liiXl)qi{X2) + qi{Xi)qi{X2)} 


(19) 


The asymmetry is significant only when xi is large and 
0^2 is small (when X 2 is small, u{x 2 ) ~ d{x 2 ) ~ u{x 2 ) ~ 
d{x 2 )). The asymmetry for u quarks dominates, and the 
fractions of d quarks and u antiquarks are sources of dilu¬ 
tion. 


D 


LHC 

AFB 


{d)?S 


D 


LHC 

AFB 


{q)?S 


d{xi)d{x2) 

u{xi)u{x2) 

u{xi)u{x2) 

u{xi)u{x2) 





( 20 ) 

( 21 ) 


Since xi = both the mass and rapidity depen¬ 

dence of Afb provides information on PDFs. 


At the LHC, the W asymmetry also provides informa¬ 
tion on the d/u ratio. The W~ /W^ ratio at the LHC can 
be written as 


/Y^\LHC 


d{xi) u{x2) ± s(xi) c{x2) 
u{xi) d{x2) ±c(xi) s{x2) 
d/u{xi) ^ d 
d/u{x2) U ^ 


( 22 ) 


Unlike the situation at the Tevatron, more precise W 
asymmetry measurements at the LHC provide informa¬ 
tion on the absolute value of ^(^i). Therefore, new mea¬ 
surements of the W charge asymmetry at the LHC (which 
have not yet been incorporated into PDF fits) can be used 
in combination with the constraints from Afb to reduce 
the PDF uncertainty in the extractions of sin^^e// and 
sin^ Ow at the LHC. 

Combining constraints from both Afb and new W 
asymmetry measurements can be done by adding the val¬ 
ues of Xwasym the comparison of the new W asym¬ 
metry data with the predicted W asymmetry for each PDF 
replica, to the x\fb values from the fits to extract sin^ 6>e// 
from the AFB{M,y) data for each PDF replica. 


4.1 Mass dependance of AFB{M,y) as a function of 
siH^Oeff and PDFs at the LHC 

The top panel of Fig. shows AFB{M,y) at the LHC at 
yds=8 TeV for six rapidity bins 0 < |^| < 0.4, 0.4 < \y\ < 
0.8, 0.8 < 1^1 < 1.2, 1.2 < 1^1 < 1.6, 1.6 < \y\ < 1.0 and 
2.0 < 1^1 < 2.4 (iY=0 to 5). These six bins have average 
1^1 values of 0.2, 0.6, 1.0, 1.4, 1.8 and 2.2. The mass bins 
are 60-70, 70-78, 78-84, 84-87, 87-89, 89-91, 91-93, 93-95, 
95-98, 98-104, 104-112 and 112-120 GeV. The horizontal 
scale for each of the six plots is the dimuon invariant mass 
for each rapidity bin expressed as 12 x zY ± iMass. 

The calculations are done with the POWHEG [17] MC 
generator. The version of POWHEG that is used does not 
include electroweak radiative corrections. Therefore, this 
version of POWHEG requires an input value of sin^ Oef / for 
the calculation of Afb 

The green bands in the bottom panel of Fig. span 
the difference between Afb{M, y) calculated with the 100 
NNPDF 3.0 (nlo) replicas and AFB{M,y) calculated with 
the default nnpdf 3.0 (nlo) PDF. 

The blue lines are the differences between AFB{M^y) 
calculated for several values of sin^ Oeff (sin^ 6>e//=0.23120 
±0.00040, ±0.00080 and ±0.00120) and AFB{M,y) for 
the nominal sin^ ^e//=0-23120. For all of the blue lines, 
Afb{M^ y) is calculated with the default nnpdf 3.0(nlo) 
PDF. 

As is the case for the Tevatron, the dependence of 
AFB{M^y) on siA^Oeff and on PDFs is different. In the 
region of the Z pole, AFB{M,y) is sensitive to the vec¬ 
tor couplings, which are functions of sin^6>e//. At higher 
and lower mass Ai?^(M, ^) is sensitive to the axial cou¬ 
pling and therefore insensitive to value of sin^ Oef /. As is 
the case for the Tevatron, the magnitude of the dilution 
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Fig. 8. Analysis of one of the 64 LHC pseudo experiments 
(6.7 M dimuon events with CMS-like detector acceptance cuts) 
with 100 PDF replicas. The pseudo data are generated by 
the POWHEG MC with the default NNPDF 3.0 (nlo) PDF and 
sin^ ^e//=0.23120. The top two panels show the extracted 
sin^ Oeff and corresponding values from fits to Afb{M, y) 
versus replica number for the 100 NNPDF 3.0 (nlo) replicas. 
The bottom panel shows the same results in the form of a 
scatter plot of Xa/ 6 values versus sin^ Oeff for one pseudo ex¬ 
periment. The number of degrees of freedom is 71 (=6 x 12 — 1). 


of Afb{M) is larger in regions where the absolute value 
of Afb{M) is large (i.e. away from the Z pole). At the 
LHC the dilution depends on both M and y. The com¬ 
bined mass and rapidity dependence of the dilution at the 
LHC provides more stringent constraints on PDFs than 
Afb{M) measurements at the Tevatron. 


4.2 MC studies with NNPDF 3.0 PDFs at the LHC 

For studies of AFB{M^y) at the LHC we simulate Drell- 
Yan dimuon data for 64 pseudo experiments for a CMS like 
detector at \/s=8 TeV. The pseudo data are generated by 
the POWHEG NLO MC generator with the default NNPDF 
3.0 (nlo) PDFs, The pseudo data are generated with an 
effective mixing angle sin^ ^e//=0.23120. 

For each pseudo experiment, we generate a sample of 
15.6 Million dimuon events with > 50 GeV, which 
corresponds to an integrated luminosity of 15.0 fb“^. This 
is similar to the ~19 fb“^ of integrated luminosity col¬ 
lected by CMS and ATLAS at 8 TeV. We apply accep¬ 
tance and transverse momentum cuts which are similar to 
a CMS-like detector. We also smear the muon energy with 


a muon momentum resolution similar to a CMS-like detec¬ 
tor. The final sample consists 6.7M reconstructed dimuon 
events. 

The 8 TeV W decay lepton asymmetry data at the 
LHC has not yet been incorporated into the most recent 
PDF fits. Therefore, in addition to Ai7^(M, ^), we also 
use the default NNPDF 3.0 (nlo) and generate pseudo data 
for the W muon decay asymmetry as a function of muon 
rapidity (for muon transverse momentum PT>25 GeV). 
This simulates the W asymmetry measurement at 8 TeV. 

In the analysis of each of the 64 pseudo experiments 
generated with the default NNPDF 3.0 (nlo) the extracted 
values of Ai?^(M, ^) for each experiment are compared to 
^) templates. The templates are generated with 
the POWHEG MC for a range of values of sin^ Oeff for each 
of the 100 NNPDF 3.0 (nlo) PDF replica. For each replica 
we extract the best fit value of sin^ ^e//, the corresponding 
statistical error and the fit XAfb- 

In addition, we calculate Xwasym which is the 
the agreement between the predictions for the W lepton 
decay asymmetry and the W lepton decay asymmetry 
pseudo data at 8 TeV for each of the 100 PDF replicas. 

Fig-i shows the results from one of the 64 pseudo 
experiments at the LHC. The top two panels show the 
extracted sin^ Oef / and corresponding XAfb values from 
fits to Afb{M^ y) versus replica number for the 100 NNPDF 
3.0 (nlo) replicas. The bottom shows the same results in 
the form of a scatter plot of XAfb values versus sin^ Oef / for 
one pseudo experiment. The number of degrees of freedom 
is 71 (=6 X 12 - 1). 

For each pseudo experiment we find the mean value 
and PDF uncertainty of sin^ Oef / from the average and 
RMS of the sin^ Oeff for the 100 PDF replicas. The average 
and RMS values are done in three ways: 

1. Using the standard average and RMS of the Oeff 
fit values. This analysis results in a standard PDF un¬ 
certainty of ±0.00051 with 100 replicas. 

2. Using the XAfb values of the fits to Afb{M, y) to form 
a weighted average and weighted RMS of the sin^ Oef / 
values. This analysis results in a PDF uncertainty of 
±0.00029 with 37 effective replicas. 

3. Using the combined x\fb^Xwasym Drell- 

Yan Afb{M^ y) pseudo data and the fits to the W lep¬ 
ton decay asymmetry pseudo data to form the weighted 
average and weighted RMS of the sin^ Oeff values. This 
analysis results in a PDF uncertainty of ±0.00026 with 
15 effective replicas. 


4.3 Studies with 1000 replicas 

As shown in Table the number of effective PDF repli¬ 
cas is reduced to 15 when we apply constraints from both 
Xj^fh RRd Xwasym' PDF Uncertainty is reduced to 
±0.00026. The uncertainty in the estimate of the PDF 
uncertainty is ±0.00005. If we start with 1000 PDF repli¬ 
cas, the number of effective PDF replicas is ~150, and 
the uncertainty in the estimate of the PDF uncertainty is 
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Fig. 9. Scatter plots of Xa ft values versus sin^ Oef / for one 
of the 64 LHC pseudo experiments. Here templates are gener¬ 
ated with 1000 replicas for (a) NNPDF 3.0(nlo) (b) CTIO(nlo), 
(c) CT14(nlo), and (d) MMHT(nlo). The number of degrees 
of freedom is 71 (=6 x 12—1). The pseudo data are gener¬ 
ated with POWHEG with the default NNPDF 3.0 (nlo) PDF and 
sin^ ^e//=0.23120. (6.7 M dimuon events with CMS-like detec¬ 
tor acceptance cuts). 


Fig. 10. The average of the results from the analyses of the 64 
LHC pseudo experiments. Each pseudo experiment is analyzed 
with 100 NNPDF3.0 templates, 1000 NNPDF3.0 templates, 
1000 CTIO templates and 1000 MHHT templates. The pseudo 
data for each experiment are generated by the POWHEG MC 
with the default NNPDF 3.0 (nlo) PDF and sin^ ^e//=0. 23120. 
(a) Analysis using the standard mean and RMS of the sin^ Oeff 
values extracted with each PDF set. (b) Analysis using the 
XAfb weighted mean and RMS of the sin^ Oeff values extracted 
with each PDF set. 


reduced to ±0.00002. Therefore, the analysis is somewhat 
more robust if we start with 1000 PDF replicas. 

Fig.^shows scatter plots of XAfb values versus sin^ Oeff 
for one oi the 64 LHC pseudo experiments. Here templates 
are generated with 1000 replicas for (a) NNPDF 3.0(nlo) 
PDF set (b) CTIO(nlo) PDF set, (c) CT14(nlo) PDF 
set and (d) MMHT (nlo) PDF set. The number of de¬ 
grees of freedom is 71 (=6 x 12 — 1). The pseudo data 
are generated with POWHEG with the default NNPDF 3.0 


(nlo) PDF and sin^ ^e//=0-23120. (6.7 M dimuon events 
with CMS-like detector acceptance cuts). 


In order to reduce the statistical error and investigate 
the PDF uncertainties, we take the average of 64 pseudo 
experiments. The statistical error in the average of the 64 
sin^ Oeff measurements is ±0.00007 (=0.00052/8). Fig. 
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shows the average of the results from the analyses of all 64 
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Table 3. Values of sin^ 6w with statistical errors and PDF 
uncertainties expected for a 15 fb“^ Drell-Yan dimuon sample 
at the LHC (at 8 TeV). The pseudo data are generated by the 
POWHEG MC generator with the default NNPDF 3.0 (nlo) PDF, 
and sin^ ^e//=0.23120. The PDF uncertainty for a standard 
analysis is compared to the PDF uncertainty for an analysis 
with both XAfb weighting and x\fb + Xwasym weighting. 


LHC CMS like 

Pseudo-Experiment 

LHC 15 fb“^ 8 TeV 

6.7M A 

reconstructed events 

input 

POWHEG 

Default 

NNPDF 3.0 (nlo) 
(261000) 

sin^ Oef f input 

0.23120 

statistical error 

A sin^ Oef f 

CTIO PDE error 

±0.00050 

± 0.00080 

Analysis replicas 

NNPDF set 

Templates 

100 

NNPDF 3.0 (nlo) 
POWHEG 

Average method 
extracted sin^ Oef / 
Standard PDE error RMS 
(uncertainty in PDE error) 

Neff = 100 

0.23121 
±0.00051 
( 0.00004) 

XAfb weighting 
extracted sin^ Ow 
x\fh weighted PDE error RMS 
(uncertainty in PDE error) 

Neff = 37 
0.23119 

±0.00029 
( 0.00003) 

XAfb+Xwsym Weighting 
extracted sin^ Ow 
Weighted PDE error RMS 
(uncertainty in PDE error) 

Neff = 15 
0.23122 
±0.00026 
( 0.00005) 


LHC pseudo experiments with templates generated with 
100 NNPDF3.0 replicas, 1000 NNPDF3.0 replicas, 1000 
CTIO replicas and 1000 MHHT replicas. The standard 
mean and RMS (=PDF uncertainty) of the sin^ Oe fj v al- 
ues extracted with each PDF set are shown in Fig. P^a). 
The XAfb weighted mean and RMS(=PDF uncertainty) 
of the sin^ Oeff values extracted with each PDF set are 
shown in Fig. H^b). 

As expected, since the pseudo data are generated with 
POWHEG with the default nnpdf 3.0 (nlo) PDF, the in¬ 
put value of sin^ ^e//=0.23120 is extracted with no bias 
when the pseudo data are analyzed using templates gener¬ 
ated with either 100 or 1000 nnpdf 3.0 (nlo) replicas. The 
PDF uncertainty is reduced from ±0.00052 to ±0.00030 
when XAfb weighted mean and RMS are used. 

The CTIO PDFs are less precise because they do not 
incorporate any LHC data. Consequently, the uncertain¬ 
ties with CTIO PDFs are larger. The CTIO PDF uncer¬ 
tainty is reduced from ±0.00078 to ±0.00036 when x\fb 
weighted mean and RMS are used. Similarly, the bias 
with CTIO is reduced from ±0.00031 to -0.00026 which 
is within the reduced PDF uncertainty. The CT14 PDFs 


Table 4. Expected statistical and weighted PDFs uncertain¬ 
ties in the measurements of sin^ Ow and with a CMS 

like detector for two samples, (a) A total of 15M reconstructed 
dilepton events (8.2 M and 6.8M e+e“) in a CMS like 

detector. This is is similar to the existing 19 fb“^ CMS data 
sample at 8 TeV. (b) 120M reconstructed events, which 

the sample expected for a CMS like detector with 200 fb“^ at 
13-14 TeV 


CMS like detector 

2016 

2017-18 


sample 

sample 

Energy 

8 TeV 

13-14 TeV 

Number of 

8.2M AlJ‘~ 

120M n+IJ'~ 

reconstructed events 

6.8M e'’'e“ 

- 

A sin^ Ow 



Statistical error 

± 0.00034 

± 0.00011 

Weighted PDE error 

± 0.00022 

± 0.00014 

(Stat±PDE) error 

± 0.00040 

± 0.00018 

^M^ndirect 

MeV 

MeV 

Statistical error 

±17 

±5 

weighted PDE error 

±11 

±7 

(Stat±PDE) error 

±20 

±9 


and MMHT PDFs incorporate LHC data in the fits. The 
PDF uncertainties with CT14 are reduced from ±0.00051 
to ±0.00034 when x^j^fb weighted mean and RMS are used. 
Similarly, the bias with CT14 is reduced from ±0.00022 
to -0.00016, which is within the reduced PDF uncertainty. 
The PDF uncertainties with MMHT are reduced from 
±0.00051 to ±0.00029 with x\fb weighted mean and RMS. 
Here, the bias with MMHT is reduced from -0.00063 to - 
0.00044, but it is still larger than the PDF uncertainty. 

As shown in Fig. [9, the analysis of the pseudo data 
illustrates that MMHT PDF set is not fully consistent 
with the NNPDF or with CT14 PDF set. A similar study 
with actual data at 8 TeV would be a first step in the 
investigation of the origin of the differences between the 
various PDF sets. 


5 Conclusion 

We show that measurements of the Drell-Yan forward- 
backward charge asymmetry (Ai7^(M, ^)) at hadron col¬ 
liders provide a new powerful tool to reduce the PDF un¬ 
certainties in the measurement of electroweak parameters. 

Table summarizes the analysis for two samples. The 
first (labeled 2016) is a sample of 8.2M and 6.8M 

e+e“ reconstructed events (with Mu > 50 GeV) corre¬ 
sponding to an integrated luminosity of 19 fb“^ for a CMS 
like detector at 8 TeV. This sample is similar to the ex¬ 
isting 19 fb“^ CMS data sample at 8 TeV. The statistical 
error in the measurement of sin^ Oef f for this sample is 
expected to be ± 0.00034, and the weighted PDF uncer¬ 
tainty is expected to be ± 0.00022. These are equivalent 
to a statistical error of ±17 MeV and a weighted PDF 
uncertainty of ±11 MeV in the indirect measurement of 
Mw • 
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With the larger number of events expected to be 

collected at 13-14 TeV, both the statistical errors and the 
weighted PDF uncertainties are expected to be smaller. 
About 120M reconstructed events (with > 50 

GeV) are expected in a CMS like detector for an integrated 
luminosity of 200 fb“^ at 13-14 TeV. For this sample (la¬ 
beled 2017-18), as shown in the second column of Table 
the expected statistical error in the indirect measurement 
of Mw is 5 MeV, and the weighted PDF uncertainty is ±7 
MeV. These expected errors are smaller than the uncer¬ 
tainties in the most recent direct measurements of Mw- 
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